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We study a trial wavefunction for an interstitial in a Wigner crystal. We 
find that the electron correlations, ignored in a conventional Hartree-Fock 
^J ■ treatment, dramatically lower the interstitial energy, especially at fillings close 

f~^ ■ to an incompressible liquid state. The correlation between the interstitial 

^: 

0^ . electron and the lattice electrons at i/ < 1/m is introduced by constructing 



O 



1 



a trial wavefunction which bears a Jastrow factor of a Laughlin state at v 



B 

i-rt ■ 1/m. For fillings close to but just below z^ = 1/m, we find that a perfect 

O ■ Wigner crystal becomes unstable against formation of such interstitials. It is 



argued that conduction due to correlated interstitials in the presence of weak 

disorder leads to the classical Hall resistivity, as seen experimentally. 
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In very strong perpendicular magnetic fields at low temperatures, it is generally believed 
that a two-dimensional electron gas should form a crystalline lattice, in the absence of dis- 
order. A growing body of experimental evidence has accumulated in the last five years that 
this so-called "Wigner crystal" |I]] (WC) may have been realized in high-quality heterojunc- 
tion systems. Evidence of this has come in the form of rf data @J§], transport experiments 
10], cyclotron resonance p, and photoluminescence (PL) experiments [P,|3|- One of the 
very interesting aspects of this system is the competition between the quantum mechanical 
correlations that lead to the fractional quantum Hall effect (FQHE), which are elegantly 
described by the Laughlin trial wavefunctions |^, and the types of wavefunctions that de- 
scribe the Wigner crystal, which is typically described by a Hartree or Hartree-Fock (HP) 
approximation as a starting point 0. In particular, transport data [Q have suggested that 
there is a very sharp transition between the PQHE and the WC as one varies the magnetic 
field in the vicinity of filling fraction v = 1/5. Since the transition occurs over an extremely 
narrow magnetic range, it is natural to ask whether there are any fiuctuations in the crystal 
state that hint at the transition if one is at filling fractions (for example) near, but below 
the transition. We address this subject in this Letter. 

Our principle motivation for this study is the expectation that density fiuctuations of the 
WC in the vicinity of a PQHE filling fraction should be strongly infiuenced by Laughlin-like 
correlations. As a paradigm of this, consider the WC at filling fractions close to, but just 
below u = 1/5. Among the basic excitations of the WC, interstitial electrons in particular 
create a local charge density that is raised closer to, or possibly above, 1/5. We expect that 



interstitial energies should be strongly renormalized by correlations in this situation [jTO 
To model such effects, we consider the following trial wavefunction for a WC with a single 
interstitial: 

N 

^ = A{Y[iu - z.)'"e-l"l'/Viyc(^i, • • • , ^tv)} (1) 

where ipwci^i, . . . , z^) is the wavefunction for a perfect WC of N electrons, with elec- 
trons localized near sites R°, Zi = Xi — iyi is an electron coordinate in complex notation, u 



is the interstitial coordinate, m is an integer chosen such that 1/m is larger than the fill- 
ing fraction u (we show below that the wavefunction no longer represents an interstitial for 
z/ > 1/m), and A is the antisymmetrization operator. We work here in units of the magnetic 
length lo = {hc/eBy^"^, where B is the applied magnetic field. We note that if one uses a 
Hartree approximation for ipwc, the case m = corresponds to an antisymmetrized Hartree 
approximation for the interstitial. Such approximate forms for wavefunctions of the WC at 
low enough fillings have been argued to be quite good, since one may show that exchange 
corrections to the WC energy at low fillings are quite small [^. 

We have evaluated the energy of the wavefunction in Eq. (|l]) for a Coulomb potential 
approximately for various fillings and choices oim. Fig. || illustrates these energies for m = 5 
and m = 7 as a function of z/. The energy e plotted here is defined as e = -E^v+i ~-^]v+i5 where 
i^^^r+i is the energy of the state with N electrons and 1 interstitial, E^^^ is the energy of a 
perfect lattice of A^ + 1 electrons. Thus, e represents the energy to deform a perfect lattice of 
A^ + 1 electrons into a lattice with A^ electrons and one interstitial, with no change in the net 
charge of the system. For simplicity, we have used a square rather than a triangular lattice 
for ipwc'i this should have no effect on our basic conclusions. In the insert of Fig. |T], the 
energy e is shown as a function of the strength of the correlation parameter m. One can see 
in the main figure that the energies are extremely low over a large range of filling fractions, 
and clearly become negative for filling fractions just below u = 1/m. (For the case m = 5, 
the energy is just slightly negative over a surprisingly large range of u. Our numerical error 
unfortunately does not allow us to conclude that the energies are actually negative except 
above z/ ^ 0.17). We have found similar results for u = 1/9. 

Several important conclusions immediately follow from these results: (1) A perfect 
Hartree- Fock WC cannot be stable for filling fractions very close to z^ = 1/m, for m < 9. Our 
calculations indicate that, at the very least, a superlattice state of correlated interstitials 
would have lower energy than a perfect HF WC at such fillings. We will argue below that 
simply including a Jastrow factor in front of an otherwise perfect HF WC will not lead to 
such a large lowering of the energy. Such a structural phase transition may be related to PL 



anomalies seen for v = 1/m, but at much lower fillings than those for which the FQHE is 
found in transport [jlT|. (2) Since the energies of the interstitial states are so very low over a 
large range of filling fractions, these are likely to contribute significantly to charge transport 
at fillings for which the crystal is stable. Furthermore, we will argue below that the Hall 
resistivity in the presence of weak disorder takes on the classical value p^y = h/ue^, as is 
found in experiment, [l^ (3) Experimental observation of a vanishing activation energy as 
the filling fraction approaches 1/5 from below could be understood in terms of the correlated 



interstitials, whose energies vanish near the FQHE filling fractions. (4) The well-known [T^ 
"dip" in pxx near u = 1/7, where to date no FQHE has been observed, might be understood 
as being due to correlated interstitials, which have particularly low energies near this filling. 
To understand why the correlated interstitial has such low energy, it is convenient to use 
a variation of the "plasma analogy" first exploited by Laughlin in the context of the FQHE 



. Dropping the antisymmetrization [0 in Eq. (||), and using a Hartree approximation 



for ipwcy we write the square modulus in the form |^p = e ^*, where, for /3 = 1, 

^ |7,|2 ^ U. _ fi>o|2 

$ = -2m;^ln(|«-z,|) + ^ + ;^^^-l^ (2) 

where R" = R°^ — iR°y are the perfect WC lattice sites in complex notation. This represents 
the energy of a classical two-dimensional Coulomb particle with charge 1, interacting with a 
neutralizing background of charge density a^ = l/27r, and a lattice of charges m with charge 
density (Xi = m/a"^, where a is the WC lattice constant. On a coarse scale, the interstitial 
is thus interacting with a system of net charge density a = (1 — r7iz/)/27r (z/ is the filling 
in the absence of the interstitial), and thus spreads out in the center of the disk to an area 
of l/o". We note that precisely at the filling u = 1/m, the interstitial will "spread out" 
uniformly over the disk, and for higher fillings the interstitial electron will have its highest 
probability at the system edge, essentially being ejected from the system. This is the reason 
the wavefunctions are only physically reasonable for fillings below u = 1/m [Q. Fig. |^ 
illustrates the probability of finding the interstitial at a particular position along a crystal 
symmetry axis for ?7i = 5 and various values of z/, which is calculated using Eq. (^ below. 



We are now in a position to understand why the correlations lower the energy of intersti- 
tial, especially close to z/ = 1/m. In essence, the Jastrow factor allows us to take advantage 
of quantum fluctuations in the WC When there are configurations in ip^c ^^^^ have large 
holes in them, the interstitial has a high probability of being located there, especially when 
V ~ l/m, for which the interstitial "density" (Fig. 0) samples a large region of the crystal. 
Conversely, when the crystal forms regions of high density, the interstitial has a low proba- 
bility of approaching that region. In this way, the interstitial electron does an excellent job 
of minimizing its energy. We note that a Jastrow factor of exponent m, involving all electron 
pairs, multiplying a perfect WC HF state would probably not lower the energy per electron 
so much as the interstitial electron energy is lowered when v < 1/m. This is because in the 
former case the electrons tend to remain closely tied to their lattice sites, and hence do not 
"sample" enough of the lattice to take great advantage of the quantum fluctuations. Thus, 
the low energy of the interstitial wavefunction is due as much to its locally high density as 
it is to the good correlations it has with the lattice electrons. 

We now consider the Hall effect due to correlated interstitials in the presence of weak 
disorder, and show this gives the classical Hall resistivity, as seen in experiment. |]T^ We 



first need to determine how disorder is likely to affect our interstitial wavefunction. The 
main source of disorder in heterostructures is thought to be caused by fluctuations in the 
neutralizing background charge density, which typically is set back from the electron gas by 
several thousand Angstroms. Recent classical calculations of the electron crystal state in 
the presence of such disorder have shown that this leads to small variations of the crystal 



density over long length scales. [T^ As seen in Fig. 1, the interstitial energetically favors 
regions of higher density. We can build correlations into the wavefunctions to optimize the 
energy by altering our correlation factor slightly, to take the form ]\i{u — Zi)"^\ with rrii 
chosen to be larger for lattice sites in low density regions, and smaller elsewhere. One may 
show that the interstitial will not be ejected from the system, provided < rrti >< 1/z/, where 
< . . . > denotes an average over lattice sites. The higher energy regions may be optimally 
avoided if one chooses < rrii >= 1/z/. 



The Hall resistivity of interstitials in such states can be estimated using a mean-field 
approximation, in which we replace rrii hj < nii > for all the lattice sites. This is justified, 
because if the interstitials are fiowing in a current, they will pass many of the lattice electrons. 
Furthermore, the correlation factors turn out to affect the lattice electrons over very long 
distances, as will be described below. The key insight is that the correlation factor {u — z.j)'^^ 
makes the interstitial appear as a magnetic flux tube to lattice electron 2, with strength 
rriiCpQ, where 0o = hc/e. Thus, on average, a single correlated interstitial flowing through 
the system should have the same effect as a moving flux tube of strength < rrii >= \jv. It 
is easy to show using Faraday's law that an electric current / flowing in the x direction in 
which each carrier has this magnetic flux will generate a force on the lattice electrons, which 
may be described by an effective voltage drop of magnitude -^l- This results in a flnite 
Hall resistivity p^y = -t-, as observed experimentally. [|17| We note that a similar mechanism 



leads to the flnite Hall resistance of the recently proposed "Hall insulator" . |]T9 

We now discuss in some detail how the energy of the interstitial state is evaluated. Using 
Hartree approximation for V^vFC in Eq. (0) and neglecting the anti-symmetrization, one has 

N 

\^[u, zi, . . . , z^)\' = e-l"l^/2 l[\u- ^,r-e-l--^°l^/^ (3) 

1=1 

where R° forms a square lattice, which would be the center of mass of each lattice electron if 
the Jastrow-type correlation factor in the wavefunction were absent. The interstitial electron 
density distribution p{u) = JYl^=id'^Zi\ip{u, zi, . . . , z^)]"^/ J (fuYl^i(fzi\ip(u, Zi, . . . , zj^j)]'^, 
shown in Fig. ^ can be directly evaluated from the above expression, since the integrations 
over the lattice electron coordinates may be computed analytically. The energy of this inter- 
stitial state -E]v+i consists of two parts: -E^v+i = -^i + E2, where Ei represents the Coulomb 
interaction between the interstitial electron and the lattice electrons, and E2 represents the 
pairwise interactions among the lattice electrons. Numerical evaluation of Ei can be carried 
out directly from Eq. (|^) |jl^ by analytical integration over the Zj's and numerical integra- 
tion over u. An exact calculation for E2, on the other hand, is highly impractical. The 
reasons are easily recognized. To obtain the interaction energy from one pair of electrons, 

6 



one needs to compute a 6- dimensional integral over the coordinates of the interstitial and 
the pair of the electrons involved. This 6-dimensional integral has to be performed for each 
of the (1/2)A^(A^ — 1) pairs for a system of N electrons. Moreover, one needs to compute E2 



for relatively large systems, since E2 turns out to depend on the system size |]T5[. We thus 
turn to an approximate evaluation of £'2- 

Our approximation scheme is based on the observation that the main effect of the Jas- 
trow factor in Eq. (|^) is to shift the center of mass of each lattice electron wavepacket 
away from R°. Beyond this, these wavepackets are also deformed from circularly sym- 
metric Gaussian distributions. Such effects can be systematically accounted for if we ex- 
press the Coulumb potential in terms of a multipole expansion. One immediately realizes 
that the effect of the center of mass shift on E2 is much more important than that of the 
wavepacket deformation. This is because the shift of the center of mass of each lattice 
electron changes the dipole and higher-order multipole moments, while the wavepacket de- 
formation changes only the quadrupole and higher- order multipole moments. Since the 
extension of the electron wavepacket is much smaller than the inter-particle distance for 
i^ ^ 1/5) a natural approximation for E2 is to neglect the wavepacket deformation, ie., to 
assume that an electron at lattice site i consists of a Gaussian wavepacket localized at the 
shifted center of mass < r^ >, where < r^ > is calculated numerically from the probability 
distribution Eq. (H) as < fj >= J (PuY[i(Pzirj\'ip\'^/ J cPuY[iCpZi\'tp\'^. Therefore, one has 
E2 — J2i<j ^(1 < "^i > — < "0' > 1); where the effective potential is taken as that of a HF 
crystal 
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The energy of the correlated interstitial state with Ei calculated from Eq. (|^) and E2 
calculated from Eq. (^ is shown in Fig. |I[ This is our primary result. Before concluding 
the description for the calculation, we remind the reader that the numerical uncertainties, 
namely, neglecting the anti-symmetrization, neglecting part of the higher-order multipole 
contribution to E2, as well as the numerical errors, are all relatively small. More important 



corrections to the precise energies of the interstitials would result from using a triangular 
lattice, and accounting for lattice relaxation in the vicinity of the interstitial. However, such 
improvements should not alter the qualitative result. 

To summarize, we have studied a trial wavefunction for an interstitial in a Wigner crystal 
and showed that the electron correlations, ignored in a conventional Hartree-Fock treatment, 
may be essential in understanding several properties of a Wigner crystal, especially at fillings 
close to an incompressible liquid state. As a result of the correlation, the state of a Wigner 
crystal with interstitials is found to have very low energy compared to the corresponding 
Hartree-Fock interstitial. In particular, we have found that for fillings close to but below 
V = 1/m, the correlated interstitial state is actually energetically favored over a perfect 
Wigner crystal. Away from such fillings, we found in the presence of weak disorder that 
thermally excited interstitials carrying a current generate a transverse electric field, leading 
to a finite Hall resistivity, pxy = h/ue'^, as observed experimentally. 

The authors would like to thank Dr. A.H. MacDonald for helpful discussion. This work 
is supported by NSF through Grant No. DMR-9202255. 
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FIGURES 



FIG. 1. 



The energy e = -Eyv+i ~ ^n^+i ^s functions of the Landau level filling factor u for m = 5 
and m = 7, where -E]v+i i^ ^^^ energy of the A^ + 1 electron state described by Eq. (^ and 
E^^i is the energy of a A^ + 1 electron Hartree-Fock Wigner crystal state. The insert shows 
the energy e at z/ = 0.19 as a function of the exponent m of the Jastrow factor, ie., as a 
function of the strength of the correlation. 

FIG. 2. 

The probability distribution for the interstitial electron in a correlated interstitial state 
described by Eq. (^) with m = 5 and z/ = 0.19,0.17, 0.1. The distance is along the y = 
axis in a Cartesian coordinate reference frame where the square lattice R°j is expressed as 
{i — 0.5, j — 0.5). The distance is in unit of the lattice constant a. 
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